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ABSTRACT 

Stellar-population analyses of today's galaxies show "downsizing", where the stars in more 
massive galaxies tend to have formed earlier and over a shorter time span. We show that 
this phenomenon is not necessarily "anti-hierarchical" but rather has its natural roots in the 
bottom-up clustering process of dark-matter haloes. While the main progenitor does indeed 
show an opposite effect, the integrated mass in all the progenitors down to a given minimum 
mass shows a robust downsizing that is qualitatively similar to what has been observed. These 
results are derived analytically from the standard extended Press Schechter (EPS) theory, and 
are confirmed by merger trees based on EPS or drawn from A^-body simulations. The down- 
sizing is valid for any minimum mass, as long as it is the same for all haloes at any given time, 
but the effect is weaker for smaller minimum mass. If efficient star formation is triggered by 
atomic cooling, then a minimum halo mass arises naturally from the minimum virial temper- 
ature for cooling, T ~ lO^K, though for such a small minimum mass the effect is weaker 
than observed. Baryonic feedback effects, which are expected to stretch the duration of star 
formation in small galaxies and shut it down in massive haloes at late epochs, are likely to 
play a subsequent role in shaping up the final downsizing behaviour. Other appearances of 
downsizing, such as the decline with time of the typical mass of star-forming galaxies, may 
not be attributed to the gravitational clustering process but rather arise from the gas processes. 
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1 INTRODUCTION 

A key issue in the study of galaxy formation is the anti-correlation 
between the stellar mass of a galaxy and the formation epoch of the 
stars in it, which is referred to in general terms as "downsizing". 
In its most pronounced form, this is simply the fact that ellipti- 
cal galaxies consist of old stellar populations and tend to be more 
massive while disc galaxies have younger stars and are typically 
less massive. However, a similar correlation between stellar mass 
and age is detected within each of the two major classes of galax- 
ies, whether they are classified morphologically as ellipticals ver- 
sus spirals or by colour as "red sequence" versus "blue sequence" 
galaxies. These trends are quite robust, e.g., they are insensitive to 
how luminosity is translated to stellar mass and colour to stellar 
age. 

A downsizing effect can actually appear in different forms 
which refer to different phenomena, involving different types of 
galaxies and different epochs in their histories. One form, which 
is the main focus of the current paper, is the fact that the star for- 
mation histories inferred from present-day galaxies using synthetic 
stellar evolution models correlate with galactic stellar mass. The 
stars in more massive galaxies tend to form at an earlier epoch and 
over a shorter time span. This pheno menon is termed "a rchaeo- 
logical downsizing" (ADS, following iThomas et aljEoO^ . Using 
observed line indices and abundance ratios, ADS has been detected 



in elliptical galaxies iThomas et ai]l2005l : iNelan et alj|2005h . and 

in a large sample of galaxies from the Sloan Digital Sky Survey 
(Heavens et al. 2004; Jimenez et al. 2005). 

The other face of downsizing is the fact that the sites of active 
star formation shift from high-mass galaxies at early times to lower 
mass systems at lower redshift. We term this phenom enon "down- 
sizing in time" (DST). It has first been detected by Co wie et alj 
( 1996), who found that the maximum rest-frame if-band luminos- 
ity of galaxies undergoing rapid star formation has been declining 
smoothly with time in the redshift range z = 0.2 — 1.7. This DST 
pheno menon has been confirmed by numerous subsequent stud- 
ies (Guzman et al. 1997; Brinchmann & Ellis 2000 ; Kodar ria et all 
■2004.: Juneau et al.. 2005; Bell et al,..2005,; Bundy et al,,200 

It is important to realize that these two forms of downsizing 
can be very different, and possibly even orthogonal to each other. 
The archaeological analysis of local galaxies highlights the forma- 
tion epoch of the majority of their stars, which at least in the case of 
ellipticals occurs at high redshifts, 2 ~ 2 — 5. In contrast, downsiz- 
ing in time refers to the specific star-formation rate (SSFR) at rela- 
tively low redshifts, 2 < 1, and therefore focuses on later phases of 
star formation, which may involve only a small fraction of the stars 
in the galaxy. Unless the stellar-mass ranking of present day galax- 
ies is the same as that of their progenitors at higher redshifts, these 
two forms of downsizing do not necessarily reflect the same phe- 
nomenon. Since hierarchical clustering in general does not preserve 
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this mass ranking, thie two forms of downsizing sliould be treated 
as two different phenomena. Indeed, as we will demonstrate below, 
the understanding of one does not imply an understanding of the 
other. 

In the standard ACDM cosmological scenario, dark-matter 
haloes are built hierarchically bottom up. This is obvious for the 
evolution of individual haloes, as they are constructed by the grad- 
ual gravitational assembly of smaller progenitor haloes that have 
collapsed and virialized earlier on. The bottom-up clustering can 
also be inferred statistically from the power spectrum of initial 
density fluctuations, which indicates that the mass distribution of 
collapsing systems is shifting in time from small to large masses. 
These hierarchical aspects of the clustering process have led to the 
misleading notion that one expects big haloes to "form" later than 
small haloes, without distinguishing between the dynamical col- 
lapse or assembly of these haloes and the formation epoch of the 
stars in them. The observed downsizing is therefore frequently re- 
ferred to in the literature as "anti-hierarchical", and thus as posing 
a severe challenge to the standard model for structure formation. 
However, when comparing the histories of different haloes, the evo- 
lution may be interpreted as bottom-up or top-down depending on 
how "formation" is defined. 

The evolution of dark matter (DM) haloes has tradition- 
ally been studied thr ough the hist ories of the main progeni tors 
jLacev & Cold 1 19931 Ei senstein Si hoeb \l99A iNu sser & Shetfl 
1999; Firmani c& Avila-Reese 2OO0I: Ivan den BoscM l2002bF 
Wechsler et al. 2002; Li et al. 2005). The main-progenitor assem- 
bly history is constructed by following back in time the most 
massive progenitor in each merger event. We term i\fmain(z) the 
main progenitor mass at redshift z. The corresponding formation 
redshift Zmain of a halo of mass Ado at z = is commonly defined 
as the time at which the main progenitor contained one half of to- 
day's mass, Af,nain(zmain) = Mq /2. Accordiug to this definition, 
more massive haloes indeed form later. The formati on redshift of 
the main progenitor has been computed bv lLacev" & Cole ( 1993) 
based on the Extended Press-Schechter (EPS) formalism, and the 
trend with mass has been confirm ed for various cosmologies (see 
for example lvan den Bosclj2002bl hereafter vdB02). This has also 
been tested using trees extracted from cosmological A'^-body sim- 
ulations (Lacev&Cole 1994; Wechsler et al. 2002). We confirm 
this behaviour below using a new analytic estimation of the full 
time evolution of A/main(2), based on the EPS formalism itself 
without the need to construct merger-tree realizations. 

However, the history of the main progenitor of a given halo 
does not represent the history of the whole population of progen- 
itors in which the stars of a present-day halo have formed. Per- 
haps more directly relevant for the stellar population at any given 
epoch is the sum over the masses of all the virialized progenitors 
in that specific tree at that time, which we term i\faii(z). If this 
summation is performed down to a zero minimum mass, we have 
by definition Alai]{z) — Ma- However, when a non-zero mini- 
mum mass Mmi-a{z) is applied, the same for all haloes, we find a 
robust archeological downsizing behaviour.^ We demonstrate this 
effect analytically based on the EPS formalism and confirm it using 
Monte-Carlo EPS merger trees as well as trees extracted from A'^- 
body simulations. We prove that this phenomenon is valid for any 
realistic power-spectrum shape and for any choice of Mmm{z), as 
long as it is the same for all haloes at a given time. We note that 



^ a similar point has been made in parallel bv lMouri & TaniguchJ i200di 
using a very different methodology. 




Figure 1. An illustration of the upsizing of Afniain versus the downsiz- 
ing of A/all in dark-halo merger trees. Compared are random trees drawn 
from the EPS probabilities for haloes of current masses Mq ~ lOMmjn 
and lOOAfnjin. The mass of the main progenitor versus the total mass 
in all the progenitors above A/njin are shown at 2 = 3. The progenitors, 
of mass M, are marked by circles of sizes and spacings proportional to 
(M/MoY^^. The values of Afmain/A/o and M^u/Mo are indicated. The 
main progenitor, along the left branch, is more massive in the less massive 
Mo, showing upsizing. The integrated mass in all the progenitors down to 
A/min is lai'ger in the massive Mq , demonstrating "downsizing". 

a similar trend has been found by iBoweil ^1991") for an Einstein- 
deSitter cosmology and a power-law power spectrum. 

The difference between Afmain and A/aii is illustrated in 
Fig.Q which compares the z = 2> progenitors above a given Afmin 
in random realizations of merger trees corresponding to current 
haloes of A/q ~ lOA/min and lOOAfmin. The downsizing be- 
haviour for A/all is apparent, while for Afniain the familiar opposite 
trend stands out (we term this trend as "upsizing"). The average 
distributions of relative masses in z = 3 progenitors, derived using 
EPS (see below) for the same two values of Afo as in Fig.Q are 
shown in Fig.|2| The upsizing of A/main is indicated by the excess 
of massive progenitors for the smaller current halo. The downsizing 
of A/all is demonstrated by the excess of the overall integral down 
to A/niin/A/o for the more massive current halo. 

A realistic and necessary condition for star formation is that 
the gas is able to cool efficiently. This is only possible if the gas 
resides in a halo whose virial temperature exceeds a critical thresh- 
old of T 10* K, above which atomic cooling becomes efficient. 
This provides a natural threshold Afniiii(z) for A/aii(z). If star for- 
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Figure 2. Upsizing of M^^^^ versus downsizing of M^n in tlie distribu- 
tion of mass in progenitors at 2 = 3. The area under eacli curve, from 
Iog(M/Mo) to 0, is the total mass in progenitors above M relative to Mq. 
The excess of mass in massive progenitors for the smaller current halo indi- 
cates upsizing of Mmain- The excess in total mass down to A/min/A/o for 
the more massive current halo demonstrates downsizing of A/aii- 



mation is of the maximum possible efficiency, namely if all the gas 
in haloes above Afmin(z) turns into stars on a free-fall time-scale, 
then the ADS in the stellar population emerges naturally from the 
ADS of A4ii(2). 

In reality, however, the star formation rate is likely to be 
slowed down by a variety of baryonic processes, especially by 
"feedback" effects. As a result, the star-formation history may 
or may not maintain the ADS seeded by Afaii(z) of the DM 
haloes. This should in principle be modeled by semi-analytic mod- 
els (SAMs) of galaxy formation, which attempt to incorporate the 
baryonic physical processes in merger trees of DM haloes. Unfor- 
tunately, early SAMs failed to reproduc e the ADS of ellipticals 
as we know it today (e.g., Kauffmann 1996; Bauah et al. 199^; 
[Kauf fmann & Chari ot 1998: Thomas 1999; Thomas & Kauffmann 
119991) . probably due to an inadequate treatment of feedback 
effects. SAMs also failed to recover the similar global trend 
obeyed by blue-sequence galaxies in color-magnitude diagrams 
ivan den Bosch 2002a; Bell et al. 2003), thus highlighting the ap- 
parent discrepancy betw een theory and observation. However , 
more recent models (e.g.. [ B ower et all2005l ; iDe Lucia et alj20oS 
ICroton et alJbOOd ICattaneo et alJl200a7 do succeed in reproduc- 
ing an ADS behaviour, largely because of an improved treatment 
of the feedback effects. The early SAMs only included super- 
nova feedback, which is efficient in slowing down star forma- 
tion preferentially i n smaller galaxies below a virial velocity of 
~ 100 km sec"^ <Dekel & Sil3ll98d) . The problem is that this 
process only causes a delay in the star formation: the gas is only 
prevented from forming stars until the halo has grown sufficiently 
massive that supernova feedback is no longer efficient. Because 
this results in relatively late star formation, even in massive galax- 
ies, the SAMs were unable to predict the correct stellar ages. 
The main success of the more modem SAMs is the inclusion of 
AGN feedback and shock heating physics, which causes a shut- 
down, rather than a delay, of sta r format ion at relati vely late times 
(e.g.. liim boim & Dekel 200 i:lBinne's| l2004; Cro ton et alJl200d 



IScannapieco et alj|2005t ICattaneo et alJlToOel; iDekel & BimboinJ 



l200d) . Although the details of AGN feedback are still poorly un- 
derstood, it has been argued that it is the main mechanism that ex- 
plains the "anti-hierarchical" nature of the relation between stellar 
mass and stellar age of galaxies. 

However, we show below that the simulated st ar formation 
histories of elliptical galaxies <De Lucia et'Zll2006h are qualita- 
tively similar to the histories predicted by Msii[z) of dark matter 
haloes. This indicates that the roots of the observed ADS can be 
found already in the natural downsizing of the dark matter haloes. 
Apparently, the complex feedback processes affecting the star for- 
mation do not change the general trend and only provide fine- 
tuning to the ADS effect. We conclude that ADS should not be re- 
garded as a surprising "anti-hierarchical" phenomenon of complex 
gas physics — it is rather the most natural, expected behaviour in 
the hierarchical clustering scenario. 

On the other hand, we find that the downsizing in time as ob- 
served at relatively low redshifts cannot be easily traced back to 
the bare properties of the dark matter merger trees. The mass dis- 
tribution of late-type efficient star formers at late times must be 
strongly affected by feedback or other gas processes and therefore 
the modeling of this aspect of downsizing should involve more re- 
alistic star formation rates. We show that only when A/min(2) is 
properly increasing with redshift, possibly mimicking the required 
baryonic effects, the star-formation rate associated with Afaii(z) 
can be forced to a qualitative agreement with the observed down- 
sizing in time. 

The paper is organized as follows. In the following introduc- 
tory section, ij2| we spell out the relevant items from the EPS for- 
malism and describe how we generate Monte-Carlo merger trees 
that serve us as a reference when needed. In ij3| we address the av- 
erage Afmain(z), dcrive au analytic approximation for it, and con- 
firm that it behaves opposite to downsizing. In Sj4|we study the aver- 
age Afaii(z), compute it analytically from the EPS formalism, and 
demonstrate that it shows a robust ADS behaviour. We also study 
the mutual correlation between the formation times associated with 
Afmain and Ms.\\. In !j5|we compute the EPS formation rate of DM 
haloes of a given mass, and compare it with star-formation histories 
in semi-analytic simulations and in observations. In Sj6|we address 
the downsizing in time of the SSFR as observed at different red- 
shifts out to 2 ~ 1. In ^ we summarize our results and discuss 
them. 

Throughout this paper we use a flat ACDM cosmology, with 
the s tandard power spect rum P{k) — kT^{k). The transfer func- 
tion iBardeen et alll98d) is 



Ml + 2.34,)^ 
^ ^ 2.34(7 

[1 + 3.89g + {IG.lqf + {5.46qf + {6.71qf 



(1) 



-1/4, 



Here q = fc/F, with k in feMpc ^ 
shape parameter JSugivamJl99l 



and r is the power spectrum 



(2) 



where fij, = 0.044 throughout the paper. Unless specifically 
stated otherwise, we use the standard cosmological parameters, 
with JIa = 0.7, r^m ~ 0.3, as = 1.0, and h = 0.7 (whenever 
we modify ilm or h we recompute F according to the above defi- 
nition). 



4 E. Neistein, F. C. van den Bosch, & A. Dekel 



2 EXTENDED PRESS-SCHECHTER THEORY 
2.1 The Formalism 

In the standard model for structure formation the initial density 
contrast (5(x) = p{-3i)/p — 1 is considered to be a Gaussian ran- 
dom field, which is therefore completely specified by the power 
spectrum P(fc). As long as (5 ^ 1 the growth of the perturbations 
is linear and (5(x,i2) = 5{'x.,ti)D{t2) / D{t\), where D{t) is the 
linear growth factor. Once 5(x) exceeds a critical threshold 5%-^^ 
the perturbation starts to collapse to form a virialized object (halo). 
In the case of spherical collapse 5%i^ ~ 1.68. In what follows we 
define So as the initial density contrast field linearly extrapolated 
to the present time. In terms of 5q, regions that have collapsed to 
form virialized objects at redshift z are then associated with those 
regions for which (5o > 5c(z) = 5''l^i^/ D[z). 

In order to assign masses to these collapsed regions, the 
Press-Schechter (PS) formalism considers the density contrast So 
smoothed with a spatial window function (filter) W{r\Rf). Here 
/?/ is a characteristic size of the filter, which is used to compute a 
halo mass M = 'yfpR'f/3, with p the mean mass density of the 
Universe and 7/ a geometrical factor that depends on the particular 
choice of filter. The ansatz of the PS formalism is that the fraction 
of mass that at redshift z is contained in haloes with masses greater 
than M is equal to two times the probability that the density con- 
trast smoothed with W{r;Rf) exceeds Sc{z). This results in the 
well known PS mass function for the comoving number density of 
haloes: 



2 _ 54z) 



da 



dM 



exp 



5c (2) 



2cr2(M) 



dM 



(3) 



jPress & SchechteJl974 . Here a'^{M) is the mass variance of the 
smoothed density field given by 



P{k) W\k;Rf 



(4) 



with W{k; Rf) the Fourier transform of W{r; Rf). 

The extended Press-Schechter (EPS) model developed by 
iBond et ai\ il99ll) . is based on the excursion set formalism. For 
each point one constructs 'trajectories' S{A'I) of the linear density 
contrast at that position as function of the smoothing mass A/. In 
what follows we adopt the notation of Lacev & Cole ( 1993, here- 
after LC93) and use the variables 5* = a'^{M) and u — 5c{z) to 
label mass and redshift, respectively. In the limit J?/ — > 00 one has 
that S — S{S) = 0, which can be considered the starting point 
of the trajectories. Increasing 5* corresponds to decreasing the filter 
mass M, and S{S) starts to wander away from zero, executing a 
random walk (if the filter is a sharp fc-space filter). The fraction of 
matter in collapsed objects in the mass interval M, M + dM at red- 
shift z is now associated with the fraction of trajectories that have 
their ^ra? upcrossing through the barrier uj — Sc {z) in the interval 
S, S + dS, which is given by 



f{S,uj) dS - 



exp 



'2S 



dS 



(5) 



72^ 

iBondetalligg T:' Bower 1991, LC93). After conversion to num- 
ber counting, this probability function yields the PS mass function 
of equation 0. Note that this approach does not suffer from the 
arbitrary factor two in the original Press & Schechter approach. 

Since for random walks the upcrossing probabilities are inde- 
pendent of the path taken (i.e., the upcrossing is a Markov process). 



the probability for a change AS in a time step Alo is simply given 
by equation J5) with S and uj replaced with AS and Au), respec- 
tively. This allows one to immediate write down the conditional 
probability that a particle in a halo of mass M2 at 22 was embed- 
ded in a halo of mass Mi at zi (with zi > Z2) as 



P(S'i,Ci;i|S2,(.^2) dSi = 
1 (cJi - CJ2) 



2^ (Si - S2)3/^ 



/(51 



exp 



S2,ti'l 
(tJl — UJ2 



iU2)dSl 
2 - 



2(Si - S2) 



dSi 



(6) 



Converting from mass weighting to number weighting, one obtains 
the average number of progenitors at zi in the mass interval Mi, 
Ml + dAIi which by redshift Z2 have merged to form a halo of 
mass M2: 

diV 

— {Mi,zi\M2,Z2) dMi = 



M2 
Ml 



P{Sl,UJl\S2,iJJ2) 



dS 



dM 



dMi 



(V) 



2.2 Constructing Merger- Trees 

The conditional mass function can be combined with Monte-Carlo 
techniques to construct merger histories (also called merger trees) 
of dark matter haloes. If one wants to construct a set of progenitor 
masses for a given parent halo mass, one needs to obey two require- 
ments. First, the number distribution of progenitor masses of many 
independent realizations needs to follow Q. Second, mass needs to 
be conserved, so that in each individual realization the sum of the 
progenitor masses is equal to the mass of the parent halo. In princi- 
ple, this requirement for mass conservation implies that the proba- 
bility for the mass of the n"^ progenitor needs to be conditional on 
the masses of the n — 1 progenitor haloes already drawn. Unfor- 
tunately, these conditional probability functions are unknown, and 
one has to resort to an approximate technique for the construction 
of merger trees. 

The most widely adopted algo rithm is the -branch tree 
method with accretion developed by ISomerville & Kolatl i 19991 
hereafter SK99). This method is more reliable than for example 
the binary-tree method of LC93. In particular, it ensures exact mass 
conservation, and yields conditional mass functions that are in good 
agreement with direct predictions from EPS theory (i.e., the method 
is self-consistent). 

The SK99 method works as follows. First a value for AS is 
drawn from the mass-weighted probability function 



/(AS, Aoj) dAS ■ 



Auj 



exp 



{Aij^ 



2AS 



dAS (8) 



AS3/2 

[cf. equation J6}]. Here Alu is a measure for the time step used 
in the merger tree, and is a free parameter (see below). The pro- 
genitor mass, Mp, corresponding to AS follows from a^(A/p) = 
a^{M) + AS. With each new progenitor it is checked whether the 
sum of the progenitor masses drawn thus far exceeds the mass of 
the parent, AI. If this is the case the progenitor is rejected and a 
new progenitor mass is drawn. Any progenitor with Mp < Afmin 
is added to the mass component A/acc that is considered to be ac- 
creted onto the parent in a smooth fashion (i.e., the formation his- 
tory of these small mass progenitors is not followed further back 
in time). Here Afmin is a free parameter that has to be chosen suf- 
ficiently small. This procedure is repeated until the total mass left, 
A/icft = M — Ma,cc — X] -^^P' is Isss than Afmin- This remaining 
mass is assigned to A/acc and one moves on to the next time step. 
As all other methods for constructing merger trees (e.g., LC93; 
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iKauffmann & White! 1993I) . the SK99 algorithm is only an approx- 
imation. In particular, it is based on the maii-weighted progeni- 
tor probability function js}, rather than on the number distribution 
Q, and mass conservation is enforced 'by hand', by rejecting pro- 
genitor masses that overflow the mass budget. Consequently, the 
number distribution of the first-drawn progenitor masses is differ- 
ent from that of the second-drawn progenitor masses, etc. Some- 
what fortunately, the sum of these distributions closely matches the 
number distribution Q of all progenitors, but only if sufficiently 
small time steps Aw are used (see SK99 and vdB02). In princi- 
ple, since the upcrossing of random walks through a boundary is 
a Markov process, the statistics of progenitor masses should be in- 
dependent of the time steps taken, indicating that the method is 
not perfectly justified. Consequently, not all statistics of the merger 
trees thus constructed are necessarily accurate, something that has 
to be kept in mind. 

In this paper we adopt a time step of 



AS 



AM 



10-3A/ 



6 + a log 



M 



(9) 



where a = 0.3, b = 0.8. As shown in SK99, this time step yields 
number distributions of progenitor masses that are in good agree- 
ment with 0. The average number of progenitors per time step 
is - 1.5 for IQ^h'^MQ ^ M lO^/i^Mf© and A/mi„ = 
IO^/i-^Mq. 
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Figure 3. Growth of the main progenitor for haloes of different present- 
day masses. The mass MmainCz) is the average at a fixed redshift z. The 
curves are normalized to match at z = 0. The halo masses, from top to 
bottom, range from 5 X 10^ to 5 X 10^^ h~^MQ equally spaced in the log. 
The symbols refer to the averages over Monte-Carlo merger trees and the 
curves represent our analytic results. The upsizing of the main progenitor is 
obvious. 



3 GROWTH OF THE MAIN PROGENITOR 

The full merger history of any individual dark matter halo is a 
complex structure containing a lot of information. It has there- 
fore been customary to define a main prof^enitor history, some- 
times termed mass accretion history I Firmani & Avila-Reese 2000 ; 
IWechsler et aL 2002, vdB02) or mass assembly history JLietajJ 
l2005h . which restricts attention to the main "trunk" of the merger 
tree. This main trunk is defined by following the branching of a 
merger tree back in time, and selecting at each branching point 
the most massive progenitor. We denote by A/main(2:) the mass of 
this main progenitor as a function of redshift z. Note that with this 
definition, the main progenitor is not necessarily the most massive 
progenitor of its generation at a given time, eventhough it never 
accretes other haloes that are more massive than itself. 



3.1 Analytical Derivation 

Using EP S merger trees and co smological A'^-body simulations, 
vdB02 and lWechsler et alJ j2002h have obtained simple fitting for- 
mulae for the main progenitor history. We show here that one 
can actually derive a useful analytical approximation for the av- 
erage Mnia.in{z), defined at each redshift as the average mass of 
Mmainiz) ovcr au cusemble of merger-trees. We derive it directly 
from the EPS formalism, without the need to construct Monte- 
Carlo merger trees. As shown in the Appendix, i\fnmin(2) obeys 
the differential equation 



Alu 



2 Mrr 



(10) 



Here S = S'(Mmain), Sq = S{Mnia.in/q), and the value of q is 
between 2 and a maximum value gmax . We show in the Appendix 
that the uncertainty in q is an intrinsic property of the EPS theory; 
different algorithms for constructing merger-trees may correspond 



to different q within the allowed range. The maximum value, Qmax, 
depends slightly on cosmology and mass. For example, gmax is be- 
tween 2.1 and 2.3 for flat cosmogonies with Qm between 0.1 and 
0.9 and halo masses between 10** and lO^'' h~'^MQ. 

Solving the differential equation for A/main we come up with 
a useful fitting formula: 



Afmain(z) 



p3 ^9 



g(32r) 

0-8 



-1^0) 



-Fa 



A/o 



.(11) 



Here g and Fq are analytic fitting functions motivated by the shape 
of the power-spectrum (see Appendix for their definition, and range 
of accuracy). For the ACDM concordance cosmology, we find that 
the standard algorithm of SK99 for constructing random merger 
trees yields an MuYa.in{z) which is well fitted by eq. <1 H with q — 
2.2. We therefore adopt this value below. Varying q between 2 and 
2.3 (the maximum range allowed) gives rise to a relatively small 
change in A^main; near Mmain = 0.5A/o this change is ~ 8%. 

Fig.|3]shows AJfmain (z) / Mq, the average, main progenitor his- 
tory for haloes of different masses today, all normalized to today's 
mass. The figure compares our analytic estimate based on equa- 
tion <1 H with the averages over histories computed from Monte- 
Carlo merger tree realizations described in i|2.2l We see that the 
analytic estimate reproduces the results from the realizations quite 
well, although there is a slight mismatch at high z. This difference 
may either reflect the allowed intrinsic uncertainty within the EPS 
formalism or it may be due to other inaccuracies in the SK99 algo- 
rithm used to construct the trees. 



3.2 Archaeological Upsizing 

We see in Fig.|3|that the average growth curve of the main progen- 
itor is shifted toward later times in more massive haloes, implying 
the opposite of downsizing, termed here as upsizing. One way to 
quantify the downsizing behaviour is via the quantity 
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Figure 4. Growth of the total mass in all the progenitors, M.^ii(z), for 
haloes of different present-day masses. The minimum progenitor mass is 
~ 10® M0, specified in equation <20l as a function of redshift. The masses, 
curves and symbols are the same as in Fig.|3] A downsizing behaviour is 
cleai'ly seen. It is more pronounced at small masses which are closer to 



1(2; I 2o,A'/o) 



di\fo 



Mo 



(12) 



Positive values of Dmain mark an archeological downsizing be- 
haviour, negative values refer to upsizing, and |Dmainl measures 
the strength of the effect. As is clear from Fig.|3| Dmain(2) < 
at all z, indicating that the main progenitor histories of dark matter 
haloes reveal upsizing. 

Is this upsizing a generic feature of Mmain(2)? To answer this, 
we write the average main progenitor mass of a halo of mass Afo a 
small time step l^uj ago as 



.(Ao.) 



Mo 



-So 



/(AS, ALj)dA5, 



(13) 



(see Appendix). Differentiating with respect to Mo while keeping 
/S.UJ fixed yields the ADS strength 



I5main(z | 20, Mo) = 

/(S,-So,Aa;) 



(14) 



^dM^^^/'^^-dM^^"^ 



Whether this is negative or not depends on the shape of S{M). For 
a self-similar power spectrum, S oc M~", we have that Dmain < 
as long as a > 0. We have also verified numerically that Dmain < 
for the standard ACDM power spectrum at all masses. While the 
above expression for Dmain is valid only for small Aci;, its sign is 
the same at all z. A more accurate expression for Dmain at any z 
can be obtained by differentiating equation <11> above. 



3.3 Assembly Time of the Main Progenitor 

Following numerous other studies (see [Q, we define the assem- 
bly redshift Zmain of a halo of mass A/o at time ujq according to 
Minain(zmain) — Mo/2. Usiug cquatiou (n) we obtain 



^(^main) — 



(15) 
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Figure 5. Formation redshifts, when the mass was one half of today's mass, 
for the main progenitor (5main. solid lines, circles) and for all the progen- 
itors (5all. dashed lines, squares) versus today's halo mass. The symbols 
and error-bars refer to an ensemble of random EPS merger-trees. The thick- 
ness of the Smain curvc refers to the allowed range obtained by valuing q 
between 2 and (jfmax in equation 1151 . The dashed line is the theoretical pre- 
diction 1231 for Sail- Zniain shows upsizing while z^n shows downsizing. 



In the case of scale-free initial conditions, were the power- spectrum 
is a pure power law, P{k) oc k", we have that S oc A/^° with 
a = (n + 3)/3. In this case the expression simplifies to 



= iUo + 



- 1) 



(16) 



We note that lLacev & Colel il993h computed a related expression 
for the average assembly redshift of the main progenitor (see the 
Appendix for more details). 

Fig. |5] shows the average assembly redshift of the main pro- 
genitor, Zmain, as a functlou of the present-day halo mass, for the 
ACDM concordance cosmology. The thickness of the curve corre- 
sponds to the allowed range of intrinsic uncertainty in q in equa- 
tion <15> , as computed in the Appendix. The ADS strength, Dmain, 
associated with the slope of 2:main(Air), does not change signifi- 
cantly with halo mass. The theoretical EPS curve shows an excel- 
lent agreement with the Zmain obtained from an ensemble of ran- 
dom EPS merger histories (circles with error bars). The error bars 
correspond to the standard deviation in Zmain over the individual 
merger-trees. 

Fig.|5]demonstrates again the upsizing behaviour of the main 
progenitor. This has been one of the reasons for interpreting the 
observed downsizing as "anti-hierarchical". 

The distribution of Zmain in our ensemble of EPS merger trees 
is plotted in Fig.|6|for three different masses. One of them is com- 
pared to the theoretical prediction by LC93, 



So 



So,co{z) — LUo)dS , 



(17) 



where / is defined in equation J5}. As discussed in the Appendix, 
the theoretical distribution agrees with the random realizations at 
low z, and any deviations are due to the limitations of the SK99 
algorithm used. 
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Figure 6. The distribution of ^main and z^n for different lialo masses. Halo 
masses are 5x lO^^^, 5x 10^^ and 5x 10^^ h'^^M^ (dotted, daslied-dotted 
and solid lines). The solid thick line is the theoretical prediction for the 
distribution of Zniaim equation 1171 . for a halo mass of 5 X 10^^ h~^MQ. 
The minimum progenitor mass is ^ 10'^ Mq, specitied in equation 1201 as 
a function of redshift. 



4 GROWTH OF ALL THE PROGENITORS 

Afniain(2) defined above only describes the mass growth history 
of the main trunk of the full merger tree. It is unlikely, however, 
that this is an honest estimator of the star formation histories of 
the associated galaxies. After all, star formation can occur in all 
progenitors that obey the necessary physical conditions, and is not 
restricted to the main progenitor. Since gas needs to cool before 
it can form stars, and since the cooling time is primarily a func- 
tion of halo mass and redshift, we assume that star formation oc- 
curs in haloes with a mass above a threshold mass, Mrain{z). This 
prompts us to define the formation history Maii (z) of a present-day 
dark matter halo as the sum of the masses of all progenitors that 
obey this condition. Supporting evidence for the possible success 
of such a model comes from the finding that the integral of A/aii(2) 
over the entire present-day halo mass function provides a useful 
backbone for understanding the observed, universal star formation 
history jHernguist & Springel 2003) . 



4.1 Analytical Derivation 

The construction of A/aii(2) for individual dark matter haloes re- 
quires a full merger tree, with a mass resolution that exceeds 
Mini-n[z). However, the formation history of a halo of mass Mo, 
averaged over many merger trees per each redshift z, can be derived 
straightforwardly from the EPS formalism. It should equal the inte- 
gral over the progenitor mass function in the range M = A/min(2) 
to Mo: 



Mall W = Mo / P{S,w\So,uo) 
V Af„,i„(z) 



AS 



AM 



AM 



(18) 



where P{S,uj\So,ujo) is defined in equation J6j. Performing the 
integral we obtain 



Mall(z) 

Mo 



erf 



V'25n.in(2) - 25*0 



(19) 



where 5111111(2) = 5[A/min(2)]. Note that Maii(2)/A/o depends 
on Mo through So, so that equation M9\ cannot be written in an 
explicit form. 

We see that for a given cosmology, the average formation his- 
tory of a halo of mass Mo is completely specified by A/inin(2). 
As a first attempt we associate Afmin with the halo mass that cor- 
responds to a virial temperature of Tvir = 10* K, the temperature 
above which atomic gas is able to cool and subsequently form stars. 
For a completely ionized, primordial gas this yields 



Mniin(2) = 1.52 X IO^'/i-'Mq ( ^ 



-1/2 



H{z) 
Ho 



,(20) 



where Avir(2) is the average overdensity of a virialized halo 
at redshift 2 relative to the critical density at that redshift 
(Brvan & Norman 1998), and H{z) is the Hubble parameter. Un- 
less specifically stated otherwise, we use this minimum threshold 
mass in what follows. 

The lines in Fig.|4|show Afaii (2) for several halo masses based 
on equation <19> . We now see that when "formation" is defined 
by Mall (2) we obtain ADS, with more massive haloes forming 
earlier. Also plotted in Fig. |4] are the results of the merger-trees 
(symbols). We see that while there is a fair, qualitative agreement 
between the histories extracted from the Monte-Carlo realizations 
and the exact EPS predictions, the level of agreement becomes pro- 
gressively worse for more massive haloes (relative to Afmin)- The 
merger trees predict an earlier formation time than what follows 
directly from EPS. A similar behaviour has been noticed by SK99 
(their Fig. 7), when comparing the empirical total mass contained 
in haloes above a minimum mass to the theoretical value. These 
deviations arise from the approximations made in the algorithm for 
constructing the Monte-Carlo merger trees (see discussion in i|2.2t . 

Figure0shows the Afaii(2) histories of individual haloes, ob- 
tained from Monte-Carlo merger tree realizations (thin lines), com- 
pared to the average formation history AJfaii(2), calculated from 
eq. 1191 (thick dashed line). The scatter is higher for lower mass 
haloes. This can be crudely understood as a Poisson noise asso- 
ciated with A'^aii, the number of progenitors above Mmin at every 
given redshift. For the massive halo. Mo = 5 X 10^^ H'^Mq, N^u 
is indeed quite large at all redshifts, leading to a small scatter. For 
the less massive halo. Mo = 5 x 10^" h^^ Mq, we have A''aii < 20 
at all redshifts, which results in a larger scatter. 



4.2 Archaeological Downsizing 

The all-progenitor histories Afaii (2) depend on the definition of the 
threshold mass Afmin(2). Here we investigate the necessary con- 
ditions for these threshold masses in order for A/aii(2) to reveal 
ADS. Similar to what was done in i)3.2l we define the "downsizing 
strength", D^u, as the derivative of Maii/Mo with respect to Mo. In 
order to study the Mo dependence, we rewrite equation J18> using 
different variables. 



Afall(t^) 



S„,i„-So 



f{AS,Lu -Loo)AAS 



(21) 



Mo JO 

where the function / is defined in equation j5j- This enables us to 
differentiate A/aii/A/o with respect to A/o, while keeping lj fixed, 
which yields 

L>aii(^ I ^0, Mo) = -/(5min - Sq,lj - Loo)-^{Mo) > . (22) 

Since / is a probability function, and AS/AAI < for all Af , we 
have that Daii is always positive. This implies that ADS occurs 
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Figure 7. Individual realizations of all-progenitor histories (thin solid 
curves) compared to the average at fixed z as calculated analytically (thick 
dashed curve). The curves are all normalized to Mq at z = 0. The upper 
and lower panels are for Mq = 5 x 10^^ and 5 X 10^'' h~^MQ respec- 
tively. Mjnin is specified in equation 1201 . 



for any choice of the threshold masses Myaini^z), and for any cos- 
mological power spectrum of fluctuations. The only assumptions 
used are (i) that the threshold is global, i.e., that i\fmiii does not de- 
pend on the specific halo mass Afo, and (ii) that the excursion-set 
trajectories are Markovian, which allows the change of variables 
leading from equation <18t to <21> . Note that the downsizing aspect 
of Maii(z) does not depend on the actual shape of /, which implies 
that ADS will occur for non-Gaussian density fluctuation fields as 
well. 

The opposite effect of upsizing could in principle occur if the 
Markovian assumption of the EPS random walks breaks down, so 
that the mass-weighted probability distribution P[S, uj\So,ujo) de- 
pends on So rather than being a function of 5 — 5o only. An ad- 
ditional requirement in this case is that the probability has a higher 
contribution from the low-mass end for larger M(, . Therefore, the 
Markovian nature of the random walks is a sufficient, but not a nec- 
essary condition for ADS to occur 



4.3 Formation Time of All progenitors 

Following the definition of assembly redshift, we define the forma- 
tion redshift of dark matter haloes, Zaii, by A/aii(2aii) = Mo/2. 
Using equation <19> we obtain 



U^ll = u{z^l{) =U0+ PVSmin - So , (23) 

where (3 = ^/2/erf(l/2) ~ 0.6745 (see also lBoweilll99lh . The 
dashed curve in Fig.|5]shows Zaii as a function of halo mass com- 
puted using equation <23> . The solid square with errorbars repre- 
sent the average and scatter as obtained from a large ensemble of 
EPS merger trees. Note that these deviate significantly from the 
direct theoretical EPS prediction, especially at large AIo- This is 
in stark contrast to the case of Zmain(A/o), where the merger tree 
results agree well with the direct theoretical predictions. This sug- 
gests that the discrepancy in 2:an(Afo) must originate in the statis- 
tics of smaller progenitors with masses < Mo/2. As shown by 
SK99, the A'^-branch tree method with accretion used for the con- 
struction of the EPS merger trees slightly overpredicts the num- 



ber of small progenitors at high redshifts. Fig. |5| shows that this 
can have a significant impact on Zaii; consequently, semi-analytical 
models for galaxy formation that are based on such EPS merger 
trees might actually overestimate the star formation rates at high 
redshifts. 

For haloes with Afo 3> Afmin we have that 2:aii > Zmain: 
typically the progenitors of a massive halo will have grown suffi- 
ciently massive to allow for star formation much before the final 
halo has assembled half its present-day mass into a single halo. 
Note that Zaii — ^main decreases with decreasing halo mass. When 
A/q ~ 2A/niin we have that Zaii ~ ^main, by definition, while 
Zaii < ^main for halocs with Afmin < A^o < 2A/fmin- Finally, 
for haloes with A/q < A/min the formation redshift Zaii is not de- 
fined. This systematic increase of ziaii — Zmain with increasing halo 
mass may have interesting implications for galaxy formation, as it 
provides a very natural means to break the self-similarity between 
haloes of different masses, and their associated galaxies. 

Although dynamical friction may delay the merging of galax- 
ies with respect to the epoch at which their host haloes merged, 
to first order we may associate Zmain with the redshift below 
which the haloes and their associated galaxies no longer experi- 
ence major mergers (i.e., below Zmain the main progenitor never 
merges with another halo of similar mass). In massive haloes, with 
Mo 3> Afmin we expect that the majority of the stars have al- 
ready formed much before these last major mergers, and this ma- 
jority of the stars will thus have experienced one or more major 
mergers since their formation. Consequently, the majority of the 
stars are most likely to reside in a spheroidal component, and the 
galaxy is an early-type with relatively old stars. Contrary, in low 
mass haloes, most of the progenitor haloes that are being accreted 
by the main progenitor at z < Zmain will have masses M < Afmin, 
and will thus not have formed stars. The gas associated with these 
progenitors can only start to form stars once they become part of 
the main progenitor: star formation and galaxy assembly occur vir- 
tually hand-in-hand, with the stars being born in-situ in what is to 
become the final galaxy at z — 0. Since the system has not under- 
gone a major merger since roughly half the stars formed, the system 
is likely to resemble a disk galaxy. 

Although this is clearly severly oversimplified, it is interest- 
ing that some of the most pronounced scaling relations of galaxies, 
namely the relations between halo mass, stellar age, and galaxy 
morphology, may well have their direct origin in the backbone of 
halo formation histories combined with a simple halo mass thresh- 
old for star formation. 

Finally, Fig.|6|shows the distribution of Zaii for haloes of dif- 
ferent masses, as obtained from our EPS merger trees. Note that the 
scatter in Zaii is smaller for more massive haloes, as expected from 
the Poisson statistics discussed in i)4.1l 



4.4 Comparison with TV-body Simulations 

While the merger trees analyzed thus far are based on the EPS for- 
malism, one can alternatively extract merger trees from cosmolog- 
ical TV-body simulations. Here the gravitational dynamics is more 
accurate, limited only by numerical resolution effects. However, it 
should be kept in mind that the identification of virialized haloes, 
and especially connecting them to construct merger trees, is a non- 
trivial enterprise involving several significant uncertainties. 

We compute Ma,n{z) from merger trees extracted from a 
ACDM cosmological TV-body simulation kindly provided by Risa 
Wechsler The simulation followed the trajectories of 256^ cold 
dark matter particles within a cubic, periodic box of comoving size 
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60 ^"'^Mpc from redshift z = 40 to the present. The particle mass 
is 1.1 X 10^ Mq, and the minim um halo mass dictated by the 
resolution is 2.2 x 10^° H'^Mq Csee lWechsler et all20o3. for de- 
tails). For the construction of M^i\{z) we impose Mmin values of 
5 X 10^° and 5 x 10^^ h~^MQ, and we compare the resulting, aver- 
age formation histories to those computed from the EPS formalism 
using the same threshold masses. The results are shown in Fig.|8| 
where symbols correspond to the formation histories extracted 
from the A'^-body simulations, while the lines show the direct, theo- 
retical predictions based on the EPS formalism (equation ll9> . Over- 
all the agreement is very satisfactory, although the A^-body simula- 
tions predict a somewhat later formation when Afmin ^ Mq. Note 
that the EPS merger trees yield formation times that are earlier 
with respect to the analytical formula (Fig.|4j- Thus, the difference 
between A'^-body simulations and EPS merger trees is larger than 
the difference between the A^'-body simulations and equation <19t . 
Despite these discrepancies, the A'^-body results clearly confirm the 
EPS prediction that A/aii of more massive haloes grows earlier. We 
therefore conclude that the ADS aspect of Maii is not an artifact 
of the EPS approximation, but is a generic property of DM merger 
trees. 

4.5 The Correlation between Formation Time and Assembly 
Time 

Since Zaii increases with increasing halo mass (ADS), while Zmain 
decreases ('upsizing'), we have that Zaii and 2:main are anti- 
correlated when considering haloes of different masses. But what 
about the relation between Zi^n and Zmain for haloes of a fixed 
mass? 

Fig.|9|shows the correlation between ^aii and Zmain for haloes 
of given masses, with Mmin = 5 x lO^'' . Results are shown for both 
the numerical simulations (solid dots) and for EPS merger trees 
(contours). For a 5 x 10^^ h~^MQ halo, the number of progenitors 
is small, and the full merger tree is not much more than the main 
trunk. As a result, the values of z^ii and Zmain are not very differ- 
ent and they exhibit a rather strong correlation. When the mass gets 
larger, the scatter in Zaii tends to zero while the scatter in Zmain re- 
mains large. Consequently, the correlation strength between Zmain 
and Zaii at fixed halo mass vanishes at large Mq. 

Th is has important im plic ations. Using a large numerical sim- 
ulation, iGao et alJ j2005h and iHarker et all j200dl have found a 
positive correlation between Zmain and the environment density: 
i.e., haloes in an overdense region assemble earlier than haloes 
of the same mass in an underdense region. If galaxy properties, 
such as stellar age, is correlated with 2:main, this means that haloes 
of a given mass host galaxies with different properties, depend- 
ing on their large scale environment. The results shown here sug- 
gest that this may be the case for relatively low mass haloes with 
Mo ~ Mmin, since these systems reveal a positive correlation 
between Zmain and Zaii. In more massive haloes, however, with 
Mo ^ Mmin, no such correlation is present, suggesting that the 
correlation between Zmain and environment will not create a simi- 
lar correlation between stellar age and environment. 

The positive correlation between Zmain and Zaii at fixed mass 
arises from their dependence on the merger tree properties at low 
redshifts. For example, assume that the merger tree for some halo is 
such that the mass al z — 0.5 is the same as at 2 = 0. In this case, 
we can use the analytical expressions for Zaii and Zmain starting at 
zo = 0.5 , and not at zq = 0. The corresponding Zmain and Zaii 
will refer to zq = 0.5, so both will be delayed by the same amount 
of time, thus establishing a positive correlation. 




Log(1+z) 

Figure 8. All-progenitor histories drawn from A^-body simulations (sym- 
bols) compared to the EPS predictions (curves). The imposed min- 
imum mass is Mmin = 5 X 10^*^ and 5 X 10^^ Mq in the 
top and bottom panels respectively. The mass bins in log mass are 
[11.62,11.78], [12.48,13.00], [13.48,14.00] where mass units are 
h~^MQ. The number of haloes within each bin is 479, 205 and 23 re- 
spectively. The EPS theoretical curves corresponding to each mass bin are 
averages over the same distribution of masses. 



When we increase the halo mass, the ratio Mmain{z) / Ms_ii{z) 
decreases (this is true for any specific redshift z). This implies that 
the fraction of mass incorporated in the main trunk is smaller, and 
as a result, there is more mass left in progenitors that belong to other 
branches. The scatter in i\faii comes from all the tree branches, 
where each branch contribute its own random behaviour. When 
Mmain(z) is Small with respect to Maii(2;) most of the contribu- 
tion to the scatter in Mail comes from branches other than the 
main. This explains why the correlation between Zaii and 2main 
gets poorer for high halo mass. 
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Figure 9. Correlations between ^main and z^\\ for random merger trees 
with the same final mass. The contours, equally spaced in the log, refer 
to the joint distribution from the EPS random merger trees. The minimum 
mass is A/min = 5 X 10^". The points are from A'^-body merger-trees, with 
the same mass bins as in Fig.lsl 



5 HALO FORMATION RATES 

We define the halo formation rate as the rate of change in Maii. 
Using equation <19> . this can be written as 



R(u)\ujq,Mq) = — 



V •S'min — So 



Mo 
exp 



(24) 



2Sn 



2So 



dui 



If we make the naive assumption that all the baryonic mass inside 
haloes with M > Afmin forms stars instantaneously, than this rate 
reflects the star formation history of galaxies that at time loo are 
located in a halo of mass Mq: these formation rates basically reflect 
the maximum possible star-formation efficiency. 

The upper panel of Fig. ^| shows the star-formation histo- 
ries of elliptical galaxies as a function of the mass of the halo in 
which these galaxies are located at z = 0, from the semi-analytical 
simulations of De Lucia et al. ( 2006, their Figure 3). Note that this 
model predicts that ellipticals in more massive haloes formed their 
stars earlier, and over a shorter period of ti me, in good, qualita- 
tive agreement with th e observational data of iThomas et alj i2005h 
and lNelan et alj i2005h . The lower panel of Fig. II 01 shows the cor- 
responding formation rates of the dark matter haloes, as defined 
by equation 1241 . Here we have adopted the same cosmological 
parameters as in De Lucia et al. (i.e., S7a = 0.75, Qm = 0.25, 
(T8 = 0.9, h — 0.73) and we have used a constant threshold mass 
of A/min ~ 1.72 X 10^" h^^ Mq, corresponding to the mass reso- 
lution of the numerical simulation used by these authors. 

The formation rates of DM haloes, as seen in Fig. 1101 reveal a 
qualitati vely similar AD S behavior as for elliptical galaxies in the 
SAM of iDe Lucia et alJ j2006) and in the observational data (e.g., 
iThomas 6^1120 05^. Although the agreement is extremely good for 
the massive haloes, the SAM predicts a significantly later forma- 
tion in lower mass haloes, indicating that the downsizing strength 
is larger in the SAM. This highlights the crudeness of our simpli- 
fied model for star formation, while assumes that stars form instan- 
taneously as soon as the halo mass exceeds A/min- The comparison 
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Figure 10. Simulated star formation rate versus EPS halo formation rate 
R for different halo masses. Top pan el: mean SSFR of ellip tical galaxies 
taken from the semi-analytic model o Lucia et all 120061 Galaxies are 
binned by their halo mass at 2: = 0. Bottom panel: maximum SFR as im- 
plied by our simplified model, namely R of equation 1241 . for different halo 
masses at 2 = 0. The curves in the two panels refer to halo masses of 10^^, 
10^^, 10^"', and 10^^ M© (solid, dashed, dotted, dashed-dotted lines re- 
spectively). Mmin is 1.72 X 10^" h~^MQ, the minimum halo mass in De 
Lucia et al. (2006). In the lower panel we add a curve for a halo mass of 
10^2 Mq with Mmin set to 10^1 Mq (thin sohd line), as an example for 
the halo formation rate when A/min is only one order of magnitude below 
the halo mass. 



with the SAM suggests that this is a fairly accurate assumption in 
massive haloes. In low mass haloes, however, the baryonic feed- 
back processes modelled in the SAM must have caused a signif- 
icant delay in the formation of the stars. Indeed, the efficiency of 
supernova feedback to cause such a delay in larger in lower mass 
haloes 1 Dekel & Silk 1986). In principle we can increase the down- 
sizing strength for the dark matter haloes by increasing A/min- For 
example, the thin solid line plots the formation rates for a halo of 
10^^ AT© but with a higher Afmin of lO" AT©. This brings the for- 
mation rates in better agreement with the specific star formation 
rates of elliptical galaxies in haloes of IO^^A/q in the SAM. Thus, 
one may mimic the delays in star formation due to supernova feed- 
back effects by an increase in the star formation threshold mass 
Afmin, eventhough we do not necessarily consider this very phys- 
ical. We conclude that the ADS in galaxies has its natural origin 
in the ADS of Afaii, while the baryonic physics associated with 
cooling, star-formation, and feedback merely causes a shifting and 
stretching of the relative formation histories. The main trend with 
halo mass, however, simply relates to the dark matter formation 
histories. 

We define the mean formation epoch of a dark matter halo as 



/g« R(uj)ujAt 
Io° Ri^)dt 



(25) 



If, for simplicity, we keep the star formation threshold mass con- 
stant, i.e., Afmin(2:) — A/min, then this reduces to 



iUR= iU0 + 



-is 

TV 



So) 



(26) 



where we have used the fact that the denominator of equation <25> 
is equal to unity. Note that this mean formation epoch is very 
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Figure 11. Effective formation epoch versus mass in EPS theory versus 
observations. The formation epoch for dark matter haloes of present mass 
Mo, based on equation |26|, is plotted for A/min = lO" and 10^^ Mq 
(solid curves). Halo masses are divided by 30 in order to roughly translate 
DM into stellar masses. The epoch for star formation as deduced from local 
ellipticals by Thomas et al. (2005) is shown (shaded area) between the two 
dashed lines which refer to ga laxies in low-density and high-density envi- 
ronments. Here h = 0.75 as in lThomas et alji20O5ft . There is a downsizing 
behaviour in both cases. 



similar to (.Jaii of equation <23> . but with /3 ~ 0.67 replaced by 
vW^i"- 0.8. 

Figure ^3 compares our analytic estimates for the mean for- 
mation epoch of DM haloes to th e star-formati on histories deduced 
from nearby elliptical galaxies bv lThomas et al. ( 2005), both for el- 
lipticals in low-density and high-density environments'^ . The solid 
lines correspond to our estimates of equation <26> for A/min = 10^ 
and 10^^ Af0, as indicated. Note that we have divided the dark mat- 
ter masses by 30 to obtain a very rough proxy for the stellar mass. 
A comparison with the data of Thomas et al. is only trully mean- 
ingful if (i) all gas in haloes with M > A/min is turned into stars 
instantaneously, and (ii) haloes host only one galaxy whose stellar 
mass is equal to Afo /30. Although neither of these is likely to be 
correct, the data and 'model' are in qualitative agreement in that the 
more massive structures have formed earlier, i.e., the model shows 
ADS. If A/min = lO'^/i"^ A/q, the DS strength is too weak, accross 
the mass range of interest, compared to the data. This indicates that 
the baryonic physics needs to delay and or suppress star formation 
relatively more in lower mass haloes. Alternatively, if A/min is sig- 
nificantly larger (~ lO^^^^A/©), the DS strength at fixed halo 
mass is stronger, and there is less need to delay or suppress star for- 
mation in order to globally match the data. However, this requires 
a yet unknown physical mechanism that can prevent star formation 
in all haloes below this mass limit. 

Yet another way to view the ADS aspect of halo formation 
histories is via the mean epoch at which a halo of mass A/q at time 
Wo has a progenitor of mass M. The number density, ANiyi), of 
progenitors with masses in the interval M to A/ + AM at time uj 
is given by equation Q. Using it to weight the averaging of lo we 
obtain 

^ The de nsity is defined as th e number of galaxies within a one degree 
radius, see lThomasetal]<2005h for details 



This resembles tjjj in equation <26> , meaning that the mean for- 
mation epoch for a given Afmin is equivalent to the mean epoch for 
progenitors to be of a given mass. A/ = A/min- The ADS behaviour 
is apparent in equation <27t from the fact that tjp increases with A/o 
(via 5*0). This implies that progenitors of a given mass appear ear- 
lier in the merger tree of a more massive present-day halo. This is 
similar to the result obtained bv iMouri & Taniguchil te006), who 
also argue that downsizing is a natural prediction of hierarchical 
formation scenarios. 



6 DOWNSIZING IN TIME 

As mentioned in ^ there is another observed downsizing effect, 
different from the archeological downsizing dealt with so far, which 
refers to the decrease with time of the characteristic mass of the 
galaxies with the highest specific star formation rates. We show 
here that, unlike the archeological downsizing, this downsizing in 
time (DST) is not in general rooted in the hierarchical clustering 
of dark matter haloes. The dark haloes show such an effect only if 
A/min is decreasing with time in a sufficiently steep pace. 

The symbols in FigfT2|show the DST obtained from the data in 
iBrinchmann & EllidOoOQy For a given SSFR we select from their 
data the stellar mass and redshift of a galaxy with that SSFR. The 
solid squares, connected by a dotted curve, plot the stellar mass of 
objects forming their stars with a SSFR of 1 Gyr~^, corresponding 
to a doubling time-scale = \ Gyr. Note that the characteristic 
stellar mass of systems forming stars at this rate is lower at lower 
redshift; this is DST. The other symbols correspond to lower SS- 
FRs of 0.1 Gyr~^ (solid dots connected by dashed curve) and 0.05 
Gyr^^ (stars connected by solid curve). Note that each of these 
curves reveals DST, and that more massive systems have lower SS- 
FRs, at each redshift. 

In order to compare this with dark matter haloes, we define 
the "current", specific formation rate of dark matter haloes as the 
rate of change of A/aii normalized to A/aii. This rate is obtained 
by setting io = ioa in the general expression for Riuj \ cjo, A/o) of 
equation j24t : 



V'5min(z)-S'(A/o) 



(28) 



For a fixed rate K — Rc, one can solve for Mc{z), the mass of 
haloes that are formed with the rate Rc'. 



S[Af,(z)] =Smin(2) 



2 /duj 
n \'dt 



(29) 



where Tc = R^^ is the corresponding time-scale. The curves with- 
out symbols in Fig 1121 show the Mc{z) relations thus obtained 
for four different time-scales Tc, as indicated. In order to allow 
for a comparison with the data, we have divided the halo masses 
by 30, as a rough proxy for stellar mass. The first thing to no- 
tice is that these 'model predictions' have almost nothing in com- 
mon with the data. First of all, all Mc{z) seem to converge to the 
same mass at low 2, independent of t^. This owes to the fact that 
/? — > 00 if A/o — > Afmin; the specific formation rate becomes 
infinite at A/min- Secondly, for high specific formation rates (low 
Tc), the Mc{z) decreases with increasing z, opposite to the DST 
observed. This simply owes to the fact that low Tc implies that 
Mc{z) ~ Mniin{z), which, according to equation <20> decreases 
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Figure 12. The mass of haloes which form at a given rate Rc at z, Mc{R = 
Rc, z) from equation I29i . The curves are marked by Tc = Rc ■ For 
high enough Tc, the mass Mc depends solely on dui/dt, while for low Tc 
it is given by Mc{z) = Mi,ji,i(2:). The connecting symbols refer to the 
observed star-formation time scale for galaxies of a given stellar mass from 
Fig. 3 of Brinchmann & Ellis i2000), for the corresponding values of Tc. 
The Mc values for the dark haloes were divided by 30 in order to allow a 
crude comparison with the stellar mass of the galaxies. 



with increasing redshift. When Tc is sufficiently high (> 10 Gyr), 
however, the dark haloes show a qualitative DST, in that Mc in- 
creases with redshift. This basically owes to the fact that the contri- 
bution from the dco/dt term in equation <29> becomes dominant 
over the term governed by the Mminiz) behavior. We conclude 
that, in general, the formation histories of dark matter haloes do 
not show a DST effect as observed for galaxies. It is clear that DST 
must be driven by baryonic processes, which must strongly decou- 
ple the star formation rates from the halo formation rates. The chal- 
lenge for the models will be to do so while maintaining a fairly tight 
coupling at high z, which, as we have shown, is required in order 
to explain the ADS. 

Another measure of DST for DM haloes is the time evolution 
of the average halo mass at which dark matter is being added to 
virialized progenitors, namely 



Me 



jR{uj\LO,M)M^dM 



rdM 



(30) 



jRii^\co,M)-^^j. 

Here 7? is given by equation <28t and dn / dM is the number density 
of haloes per comoving volume (e.g., from Sheth & Tormen 2002) . 

In Fig. 1131 we plot Ala for several different growth rates of 
Mniin{z), all normalized to coincide with our standard value of 
Mmin at z = 0. Similar to Mc defined above, the characteristic 
mass Mr is decreasing with redshift when we set Mmin to cor- 
respond to a constant virial temperature T^ir ~ 10^ K (see equa- 
tion 1201 ). or when A/min is constant in time. Only when Mmin is 
increasing with redshift roughly as 1 + 2 or faster does the mass 
Mh show a DST behaviour. 



7 CONCLUSIONS 

We have introduced a new quantity to quantify the growth of a 
dark matter halo merger tree, AIaii{z), the sum of the masses of 



Figure 13. Average halo mass Mr weighted by halo formation rate i? as a 
function of redshift for different z-dependences of M^[^{z). Downsizing 
in time is seen only when M,nin(2:) is increasing linearly with (1 + z) or 
faster. 



all the virialized progenitors at redshift z down to a minimum halo 
mass Afmin(z). We have shown, using EPS theory, that this quan- 
tity reveals an "archeological downsizing" behavior in that Ma_u{z) 
of more massive haloes grew earlier and on shorter time-scales. 
This behaviour is present for any choice of non-zero Mmin(z) and 
any cosmology. The only two conditions are (a) that the threshold 
mass Mmin(z) is independent of the mass Afo of the present-day 
halo, and (b) that the progenitor mass function, P{Mi, zi\Mo, zo) 
(equation. (6)) either depends on S{AIo) — ^(A^i) alone (i.e., the 
trajectories 5{S) are Markovian), or is such that the fraction of 
mass in progenitors below Afmin decreases with increasing A/q. 
The fact that a similar archeological downsizing effect is revealed 
by EPS merger trees and in N-body simulations indicates that these 
conditions are at least approximately valid. One should note that 
the first condition, although quite robust, might be violated in cer- 
tain circumstances. For example, today's halo mass Mq could be 
interpreted at high z as reflecting the local environment density, 
and if the threshold mass is somehow affected by its environment 
this could introduce a dependence of Mminiz) on Mo- 

Using the EPS formalism, we have analytically formulated 
the virial mass growth curve AIa_u{z), the corresponding forma- 
tion redshift Zaii, and the formation rate. The latter is found to be 
qualitatively similar to the formation rate of stars in elliptical galax- 
ies, indicating that the observed archeological downsizing in these 
systems has its roots in the formation histories of the dark matter 
haloes. However, Afaii(2:) is only a good tracer of the star forma- 
tion histories of galaxies if all the gas in haloes with M > Mmin 
forms stars instantaneously. In reality, this will not be the case, as 
cooling and various feedback processes can delay and/or prevent 
the formation of stars, even in haloes with M 3> Mmin- What is 
clear from our study, however, is that the halo mass dependence of 
these baryonic processes has to be such that it does not undo the 
mass dependence already encoded in A/aii(z). 

We have also studied the more common halo assembly histo- 
ries, defined as the mass growth histories, Mmain{z), of the main 
progenitor of the merger tree. We have developed an analytical ap- 
proximation for it based on EPS theory, and confirmed the known 
"upsizing" behaviour of this assembly history. We have shown that 
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it depends in principle on the shape of the power-spectrum, but it is 
valid for all power-law spectra as well as for the CDM power spec- 
trum. The formation times Zmain and Zaii, for a sample of equal- 
mass haloes, were found to be correlated in a way that can be un- 
derstood in terms of the mass growth at low redshifts. 

The downsizing in time, namely the decline with time of the 
mass of star-forming galaxies, cannot be easily traced back to the 
properties of the dark matter halo merger trees. With our idealized 
recipe of rapid star formation in virialized haloes above i\fmin(2), 
downsizing in time can be reproduced only if Mmin(z) is rapidly 
increasing with z. Otherwise, this kind of downsizing is most likely 
a result of feedback effects on star formation, which requires a more 
sophisticated modeling of the baryonic processes. The lesson is 
that the different faces of "downsizing" reflect different phenom- 
ena, one naturally rooted in the hierarchical dark matter clustering 
process and the other determined by non-trivial baryonic processes, 
which are yet to be properly modeled. 
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APPENDIX: ANALYTICAL FORMULAE FOR Mmain 

We use M-mai-aiz) to denote the main progenitor mass at redshift 
z, averaged (at fixed z) over many individual merger trees for the 
same parent mass Mo . Here we use the EPS formalism to derive an 
analytical estimate for Minain(2). 

Basic Equation 

Let's start with a halo of mass at time ujo, and take a small 
time-step, Au, back in time. At the time loq + Ao; we want to com- 
pute the average mass of the main progenitor. This requires the full 
probability distribution, P(Afmain|Mo, Atj), that a halo of mass 
Mo at u)o has a main progenitor of mass Mmain at time u)o + Acj. 
For A/main ^ Mo/2, one has that P(Mniain|Afo, Aoj) is equal 
to the total progenitor distribution dN /dM given by equation 0, 
simply because any progenitor whose mass exceeds Mo/2 must be 
the main progenitor. For Afmain < Mo /2, however, the only valid 
condition is that P dA'^/di\/, which is not sufficient to predict 

P(Mn,am|Mo, At^). 

As a first, naive approximation we assume that 
P(Armain|Afo, Atj) = for Afmain < Afo/2, SO that the 
main progenitor always has a mass A/main ^ Mo/2. Using 
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this approximation, tiie average mass of the main progenitor, 
Mmain(Aaj), Can be written as: 



Mmain(A(^) = / P(M|Mo, Atj)MdM. 

J Mo/2 

Using the definition of AN /AM this reduces to 



Mmam(Aaj) = Mq 



erf 



\/2S2 ~ 250 



(31) 



(32) 



with S2 = S(Mo/2) and So = S{Mo). 

We assume that Mo is just the main progenitor of the previous 
time-step"^. The rate of change, dMmain/dci;, can then be computed 
as 



AMn 



Alu 



lim 

Aui^O 



M„ 



- Mo 



(33) 



= —Mo Um - — erf 



Auj 



, V2S2 — 2So 
Using that erf (x) 2:>i/^/n when x ^ this yields 

am 



Auj 



(34) 



In the case of scale free initial conditions, were the power- 
spectrum is a pure power law (S oc M~") we can solve for Mmain 
analytically: 



(35) 



where = a [27r5(A// = 1) (2° - 1)]"^/^ 

The above derivation is based on the assumption that the main 
progenitor always has a mass A/main ^ Mq /2 in the limit Auj 
0. However, as we show below, when Auj the probability that 
Afniain < Mo/2 decrease like Auj. Consequently, this will give a 
non-negligible effect for sufficiently large uj. 



Towards better accuracy 

The dot-dashed line in Fie. 1 141 shows the distribution of the main 
progenitor masses of a halo of mass Afo ~ 10^^ Af© in a single 
time step Aui = 0.1, obtained from 10.000 realizations based on 
the SK99 algorithm. One can clearly see that P(A/maiii|A/o, Auj) 
has a non-negligible tail for A/main < Mo/2. The following anal- 
ysis aims to find the solution for A?niain taking this low-mass tail 
into account. 

The correct shape of P(A<fmain| Afo, Auj) can be constrained 
by the following conditions: 

(i) The integral of P(Afniain|A^o, Ao;) over all masses should 
equal unity, for all time-steps Auj. 

(ii) P(A/n.ain|Afo, Acj) = diV/dA/, for Afniain > Mo/2 

(equation Q), and P(A/main| Afo, Auj) < AN/AM, for Afmain < 
Afo/2 

(iii) P(A^niaiii| Afo, Alj) should not depend on the time-step 
subdivisions. This can be written as: 



P(A/"main|A/"0, At^l + AuJi) = 

J P(A/ilAfo, Aa;i)P(Af^ain|Afi, Aa;2)dAfi . 



(36) 



^ Equation 1321 becomes linear in Auj for small enough Auj, and this 
gives: Mmain(A<.i^i + A<.J2|Afo) = Afmain (Aa;2| Afmain (Aa;i| Mo)) 




M/M 







Figure 14. Average number of progenitors of mass Af one time-step Aui = 
0.1 before the present for halo of mass Afg = 10^^ h~^MQ. The solid 
curve is the theoretical prediction. The dot-dashed lines indicates that mass 
distribution of the main progenitor (defined as the most massive progenitor), 
as obtained from 10.000 Monte-Cailo EPS realizations based on the SK99 
algorithm. Note that the probability distribution for the mass of the main 
progenitor, equals dN/dM down to Mo/2, but does not vanish down to 
Af ~ 0.25J\fo. Finally, the dashed line indicates the mass distribution of 
the least massive progenitors obtained in these 10.000 realizations. 



In what follows we estimate the limits on P(Afmain| Afo, Auj) us- 
ing conditions (i) and (ii), and show that they give a narrow range 
for A/main (z). Condition (iii) does not force the solution to be 
unique, hence it enables a set of solutions, each of them is valid 
within the EPS formalism. Because condition (iii) is more difficult 
to implement, we do not compute its effect on Afmain, and assume 
that it will not significantly affect the range of solutions. 

The first condition on P(Afmain| Afo, Aoj) is that its in- 
tegral equals unity. We define ntaii as the integral over 

P(A/main|Afo, Auj) frOm Afmain = tO Afmain = Afo/2: 



ntaii — 1 



Afo 



S2 



Auj 



Af ASi-s 



exp 



Auj^ 



2AS 



AS, 



(37) 



where AS = 5" - = S{M) - S{Mo). 

We can estimate the possible effect any tail will have on 
Afmain, by Computing the effect of the most extreme tails possible. 
The first extreme is to concentrate all the tail in a small range near 
Af = 0. In this case, the integral of AfmainP(Afmain| Afo, Aw) 
over the range Afmain < Afo/2 will be zero. As a result, 
the Afmain (2) that corresponds to this extreme is given by equa- 
tion <32t . The second extreme is that all the tail is concentrated 
near Afo/2, that is, P(Afmain| Afo, Alj) has its maximum values 
(— AN /AM) down to a lower mass limit, Af/ (/max, which is set by 
the requirement that 



Mo 



Wtail = 



1 Auj 



S2 



Af ASi-5 "^""P 



Auj^ 



2AS 



AS, 



(38) 



where Sq+ = S(Af /qtmax). If we focus our attention on small time- 
steps Auj, then we can use that dP(Afmain| Afo, Atj)/dAfmain — 
near Afo/2 to approximately write that 



'9+ 



5*2 + 



TT {S2 

2 



So)' 



Auj 



(39) 
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This enables us to use a simple equation for 5,+ , combined with 
the definition of Jitaii in equation <37t . 

What remains is to find an appropriate expression for ntaii 
which is valid in the limit of small time-steps Aoj. We therefore 
split the integral in equation into two parts. The first one (m) 
is for the range < /S.S < AS^, where Alv < A^e < 1 and we 
can make the approximation AI ~ Mo: 



Alo 



= 1 - erf 



So+AS, _ 

2^ J So 

Alu 



exp 



-Au/ 
2AS 



dS 



(40) 



V2AS7 



The second range (712) is for AS > ASe where the approximation 
exp [~Auj^/{2AS)] ~ 1 is valid: 



"-2 



Mo 



S2 



1 Alo 
M ASi- 



Combining equations <37t . <40> and <41> then yields : 



Atj as: 



Mo 
2 



S2 



5o + AS< 



M ASi-s 



(41) 



(42) 



Finally we take the limit ASe ^0, and obtain 



"tail 



Alu 



^0 M - Mo dAS 
M AS^ 



+ 



\/ S2 — So 



,(43) 



Substitution in <39> then yields 



M ~ Mo dAS 
M ASi-s 



(44) 



independent of Acj. For the standard ACDM cosmology this yields 
2.1 < gmax < 2.3 for 0.1 < Qn, ^ 0.9 and lO'^/i'^Mo Mo ^ 
10^^h~^ Mq. This implies that although there is a negligible prob- 
ability that the main progenitor has a mass A/o/2.3 < i\fmain < 
Mo/2 when Au} — > 0, this probability behaves like Alo and it 
cannot be neglected. We can take this into account by rewriting 
equation <34> as 



dM„ 



dLO 



2 Mrr 



^Vs, 



(45) 



with Sq = S{Mo/q) and 2 ^ g < 2.3. In principle, any value 
of q in the range above is allowed. In particular, merger trees con- 
structed using different algorithms may have different values of q 
in the above range, as long as the algorithms adopt a sufficiently 
small time-step Alo. 

In Fig. 1151 we show A?main for a halo of mass 5 x 
lO^'^ H^^Mq. The triangles show the results obtained from many 
independent EPS merger trees, constructed using the SK99 algo- 
rithm. Note that the averaging is done over Mmain(z) at fixed 
z, which is the same as done in the analytical estimates (equa- 
tion 1311 ). The dashed, solid and dotted lines correspond to the ana- 
lytical prediction of equation <45> for q = 2, q — 2.2 and q — 2.5, 
respectively. The curve for g — 2.2 is in excellent agreement with 
the EPS merger trees. Note that this value for q is within the ex- 
pected range. In Fig.|3|we show that our analytical formula with 
g = 2.2 also accurately fits the merger-tree results for other values 
of Mo. 
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Figure 15. The mass of the main progenitor as computed in different ways 
for AIq = 5 X 10^^ h~^MQ. Big circles are from vdB02. Squares are the 
average of EPS merger trees, averaged over z at n fixed AI. Triangles are 
from the same merger trees but averaged over AI at a fixed z. Dashed lines 
show the theoretical limits of the analytic foimula (q = 2 and q = 2.5). 
The solid line is the analytic formula with g = 2.2. 



A Universal Fitting Function 

Equation <45> can be solved to obtain a direct, analytical formula 
for Mmain(z). We use the fitting function for S(M) given in 
vdB02: 



S{A1) 



cF 



1/3 



M 



1/3 



2 



g2(32r) 



where c = 3.804 x 10"'', and g{x) is an analytical function: 



g{x) = 64.087 [1 + 1.074j;" 



(46) 



(47) 



-LSSlx" ** + 0.954x°-^ - 0.185x°-'^ 



In terms of the new set of variables 



M = M- 



o(32r) 



(48) 



equation <45> does not depend on cosmology, and can be easily 
solved to give 

M = F~\tlj -U0 + FqiMo)] , (49) 
where 



F,(M) 



M Jg^{M'/q)-r{M') 



dM' 



JO M' 
Finally we write the solution in the original variables: 



Mn 



( \ ^rn rp — l 



ff(32r) 



{lo - loo) + Fq -TT-Mo 



erg ■ ■ ' 

The analytical fitting function for Fq(u) with q — 2.2 is 

Fq{u) = -6.92 X 10"^ In" ?i + 5.0 x lO"'^ In^ u+ 
+8.64 X lO'^ln^u- 12.66 In u+ 110.8, 



(50) 



.(51) 



(52) 



which is accurate to better than one percent over the range 1.6 x 

10* < AIT'^n-^ < 1.6 X 10^^ ft-^Mo 
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An Alternative Method 

We now present an alternative method to compute Mmain(2), 
which is based on the method originally introduced by LC93. 
The LC93 argument is as follows: at a specific redshift z, one 
can compute the probability for having a progenitor with a mass 
larger than Mo/2. This probability equals the probability that a 
tree will have its z-mai-n greater than z (where Zmain is defined by: 
Afniain(zmain) = Mo/2). Although LC93 claim their formula is 
only an approximation, we have not found any gap in their argu- 
ment, so we think this should be an accurate prediction. 

Lets define Q{uji\Mo,uJo) as the probability that a halo with 
mass Mo at time tJo will have its merger-tree obey Lo^Zm^i-n) = 
uj\. The probability for having a progenitor with mass bigger than 
Mo/ 2 then equals : 



^2 Mo 



So 



M 



f{S — So,(jJ - ijJo)dS ■- 



/: 



Q{LJi\Mo,uJo)duJi,(53) 



where / is defined in equation So = S'(Mo), and 5*2 = 
S(Mo/2). The distribution Q{u:>\AIo,wo) is obtained by differenti- 
ating the above equation with respect to lu and multiplying it by - 1 . 
In order to compute the mean formation time, we need to average 
Lo over the probability distribution Q: 



LjQ{Lu\AIo,ujo)duj ■ 



(54) 
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Here tJmain,i is obtained by averaging over all cj (time) possible 
for getting a mass Mo /2. This is different from the method used 
above, where Umain was computed by averaging the main progeni- 
tor masses at a fixed time. The two methods should give slightly 
different results, even if both are accurate. This is illustrated in 
Fie. 1151 where the triangles indicate the average, main progenitor 
history obtained by averaging over Mmain at fixed z, while the 
squares show the results obtained when averaging over z at fixed 

Mm am . 

So far we have repeated the analysis in LC93. Now, instead 
of computing the derivative of equation <53> . we simplify equa- 
tion <54t by a simple integration by parts: 
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/(AS, Aa;)dA5'. 



The left part is just ljo ■ We can switch the integrals on the right 
hand side of the equation, and compute the integral over w first. 
Finally we have: 

Mo /■^^-^« dAS 

Wmain,l ^ UJo = ,_ / , „ ?== ■ (5D) 



3 M{So+ AS)VAS 

For a self-similar power spectrum with S oc M~" the integral 
can be done analytically: 



Wmain.l — LOo = \ — {S2 — So) 2-F] 



1 3 



1-2" 



(57) 
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Figure 16. Mma.h-i,i(^) for several halo masses. Symbols are from the 
merger trees of vdB02. The halo masses range from 5 X 10^ to 5 X 
lQi^^h~^ Mq, spaced by a decade. Smoothed lines are the results of the 
analytic foiTnula, Eg. 1561 replacing 52 with 5(Mniain l)- The equation is 
valid only for A/main,! > Mo/2. 



where 2F1 is the Gauss Hypergeometric function. We can see that 
(ij^main.l — 1^0 )/VS2 -5o has the same value for all masses, and 
thus it is the natural variable to choose (as was done in LC93). 
On the other hand, we showed in equation <16> that (tDmain — 
<jJo)/(\/S2 — \/So) is a also constant, when the averaging of the 
trees is made along the mass axis. 

The analysis above is still valid, if we replace 5*2 with 
^(Mmain,!), SO that We cau easily generalize this result to obtain 
Mmain, 1 in the range Afmain,i > Mo /2. For masses below Mo/2, 
however, we cannot compute Mmain, i since this requires the proba- 
bility for getting a main progenitor with mass lower than Mo/2. As 
discussed above, this part of the probability function is unknown. 

In Fig. 1161 we compare results from EPS merger tree re- 
alization (vdB02, plotted as symbols) to the analytical formula 
(smoothed lines). There are some deviation between the two, pre- 
sumably because the averaging is done over a large range in red- 
shift, where the SK99 algorithm may have slight inaccuracies. This 
effect can be seen in Fig.|6|and in vdB02 (Fig. 4): the distribution 
of formation times is only accurate for low values of ci;main. Our 
previous method for computing cDmain was not affected by this in- 
accuracy because it was derived using the limit of small time-steps 
behaviour. This, combined with the fact that this method can only 
be used to compute MmsA-niz) down to Mo/2, and the fact that 
equation J56> cannot be generalized easily since 5*0 is buried inside 
the integrand, clearly favors the method discussed at the beginning 
of this appendix over the one discussed here. 



* One should take care when assigning the integral lower limit. The in- 
tegral over A5 should first be computed, similar to the analysis done for 
equation 1371 . 



